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Abstract 

By solving the Bogoliubov - de Gennes equation self-consistently in the presence 
of a non-equilibrium quasi-particle distribution, we compute the current-voltage 
characteristic of a phase coherent superconducting island with a tunnel barrier at one 
end. The results show significant structure, arising from the competition between 
scattering processes at the boundaries of the island and modification of the order 
parameter by quasi-particles and superflow. This structure is not present in non 
self-consistent descriptions of normal-superconducting nano-structures. 
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It is now well established that coherent Andreev scattering provides the key 
to understanding transport in mesoscopic superconductors and normal-superconducting 
(N-S) interfaces. For example zero bias anomaliesfl] can be understood through a de- 
scription based on multiple Andreev scattering in N-S tunnel junctions[2-5], while phase 
periodic conductances in N-S-N nano-structures[6] are understandable through theories 
of coherent transport which neglect inelastic scattering [7- 14]. All of the above theoret- 
ical descriptions are based on non self-consistent solutions of the Bogoliubov de-Gennes 
equation or corresponding quasi-classical equations and are not capable of describing the 
modification of a superconducting order parameter by a transport current. Effects of this 
kind are observable as super-gap structure in the differential conductance of N-S tunnel 
junctions and point contacts[15,16], but to date, there exists no quantitative theoretical 
framework for their understanding. For structures smaller than the inelastic phase break- 
ing length Ifi, such features cannot be ascribed to quasi-particle "heating ", because the 
energy of the electron is preserved during its passage through the sample. Instead, any 
modification is a hot electron effect and requires a description which takes into account 
the non-equilibrium distribution of electrons within the sample. 

A theoretical description which encompasses both non-equilibrium effects of 
this kind and phenomena associated with phase coherent transport does not currently 
exist. The aim of this Paper is to provide the first self-consistent description of phase 
coherent transport in a N-S-N structure, based on exact solution of the BdG equation. 
Motivated by the success of the Blonder, Tinkham and Klapwijk[17] (BTK) calculation 
for the current-voltage (I-V) characteristic of a N-S interface with a delta function scat- 
terer, we examine the simplest possible generalisation which is capable of highlighting the 
new physics which emerges from a self-consistent description. The system of interest is 
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a mesoscopic scattering region containing a superconducting island connected to perfect 
normal leads, which are in turn connected to external reservoirs at chemical potentials 
/ii and /i 2 . The system length L is assumed to be smaller than a quasi-particle phase 
breaking length and therefore a description, which incorporates quasi-particle phase co- 
herence throughout the system is appropriate. The main question of interest is whether 
or not such a description yields observable super-gap structure, which is absent from a 
BTK desciption, thereby obviating the need to introduce ad hoc heating effects. 

Before discussing quantitative results, it is useful to identify characteristic 
voltages which are missing from the BTK description [17]. For convenience, consider the 
zero temperature limit and choose /ii > /12. In this case electrons (holes) are incident on 
the island from the left (right) reservoir over an energy interval < E < fii — // (0 < 
E < ji — /i 2 ), where \i is the self-consistently determined condensate chemical potential. 
The associated current will both suppress the magnitude of the order parameter and 
generate a phase gradient. For a homogeneous superconductor with an order parameter 
A(x) = A exp[iv s x], the energy gap for excitations parallel (anti-parallel) to the phase 
gradient v s is /i_ = Ao — Pfv s (/-*+ = Aq+pfv s ), where pp is the Fermi momentum. For a 
long enough island, where the order parameter at the centre of the island is approximated 
by the above form, excitations incident on the superconductor with energies less than 
these values will be reflected and therefore in addition to the voltage A /e, one might 
naively expect the I-V characteristic to show some feature when the reservoir potentials 
satisfy /ii — /i — /i± or /i — /i 2 — H±- In addition one might expect features to occur for 
those values of /ii — /12 at which A — Pfv s = and at which the self-consistent value of 
A(x) vanishes everywhere. 

To obtain a self-consistent description, we solve the Bogoliubov - de Gennes 
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equation 

^ ^ / u n (x)\ ^ / u n (x) 
\v n {x) ) \v n (x) 

with a Hamiltonian 

/[-(fc 2 /2m)^ + «(*)-/*] A(s) 
H(x) = (2) 

V A*(x) -[-(& 2 /2m)^ + u(x)-/i] ' 

where /x is the condensate chemical potential, is the normal scattering potential 

and A(x) the superconducting order parameter, defined self-consistently by the following 
equation 



A(x) = V(x) 



S ' U n( X ) V n( X ) - ' U *n( X ) V n(x)((lL( X hn*(x))) ] (3) 

E n >0 En>0 



In this expression, primes on the sums indicate that only terms with E n less than some 
cut-off E c are to be included, due to the fact that the electron-electron interaction is only 
attractive over a small range of energies near the Fermi surface, 7^. creates a Bogoliubov 
quasi-particle and double angular brackets indicate a trace over the density matrix of 
the system. In what follows, the pairing potential V(x) is chosen to equal a constant 
for < x < L and to vanish outside this interval. The normal scattering potential is 
chosen to be u(x)//io = (2Z/kF)S(x), where /io is the condensate chemical potential in 
the absence of an applied voltage and k F = (2mfx /h 2 ) 1 ^ 2 . For a given choice of L, Z, E C) V 
and reservoir potentials, both the magnitude and phase of A(x) will be computed at all 
points in space, along with the condensate chemical potential /1. 

Since we are interested in an open system, equation (3) involves sums over 
all incoming scattering states, integrated over all E < E c . At zero temperature, for 
the case /ii > /i > /i 2 , quasi-particle states corresponding to incoming electrons (holes) 



are incident from reservoir 1 (2) over energy intervals fi± — fi (// — 112)- Assuming these 
intervals are less than the cut-off E c and if a scattering state of energy E corresponding 
to an incident quasi-particle of type a from reservoir i has a particle (hole) amplitude 
u ia (x,E) (v ia (x, E)), then equation (3) reduces to 

A(x) = V(x)j2l [ E \(ul(x,E)v l 4x,E)) + ((u* + (x,E)v l+ (x,E)))dE 

- V(x) / (u*(x,E)v 1+ (x,E))dE 

Jo 

- V(x) / (u* 2 _(x, E)v 2 _(x, E))dE (4) 

Jo 

To calculate scattering solutions in the region occupied by the island, we 
start from an initial guess for A(x) and /i and divide the interval < x < L into a large 
number of small cells of size << hp 1 , within which A(x) and u(x) are assumed constant. 
If T(xq) is the matrix obtained by producting together transfer matrices associated with 
all cells in the interval < x < x and then as outlined in appendix 1 of reference [18], the 
scattering matrix S of the island can be obtained from the transfer matrix T(L). Within 
the external leads, the most general eigenstate of H belonging to eigen-energy E is a linear 
superposition of plane waves. For a given incoming plane wave, a knowledge of S yields 
the plane wave amplitudes on the left side of the island, which can be combined with 
T(x ) to yield the wavefunction at x = x$. Given these solutions, A(x) is re-evaluated 
using equation (4) and a new choice for /i is obtained by insisting that the currents ji and 
j'2 in the leads attached to reservoirs 1 and 2 are equal. This process is repeated until the 
root mean square difference between successive order parameters is less than 1% of the 
magnitude of A(L/2). 
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In what follows all results are for an island of length Lk F = 750, a cut- 
off of E c = 0.085/io and a pairing potential of magnitude V = 0.28/io- For an infinite 
homogeneous superconductor with no current flowing and a density of states n(0), BCS 
theory predicts a bulk order parameter of magnitude A = E c /smh{l/[n(0)V]}, which 
for this choice of parameters yields A ~ 0.005/i . As an example of the results obtained, 
for an island with no barrier (ie Z — 0), figure 1 shows self-consistent results for the 
magnitude |A(x)| and phase (f>(x) of the order parameter, for various applied reservoirs 
potential differences. As expected, |A(x)| reaches a maximum value at x = L/2 and is 
suppressed at the ends of the island, on a length scale £ = kp 1 fi/\A(L/2)\, whereas the 
corresponding phase gradient v s (x) = V(f>(x) is almost a constant. Furthermore the zero 
voltage value of A (L/2), (denoted A in what follows) agrees with the BCS prediction. 
In what follows, we denote v s (L/2) and |A(L/2)| by v s and A s respectively. 

By repeating these calculations for a range of reservoir potentials and barrier 
strengths, one obtains I-V curves, whose derivative yields the differential conductance 
shown in figure 2. Clearly these curves exhibit structure which is not present in a non 
self-consistent description [17]. To identify the underlying physical processes, figure 3 
shows self-consistently determined values of A s and the various characteristic voltages 
identified above, plotted against the reservoir potential difference. For each value of Z, 
the upper and lower dashed lines of figure 3 show results for ji + and /i_ respectively, 
while the thin solid line shows A s . The upper and lower thick solid lines show values of 
Hi — fj, and (j, — /i 2 respectively. For Z — 0, where /i = (/ii — /U 2 )/2, the latter are equal. 
More generally, for the range of voltages studied, one finds /i — /i2 + «(/ii — At 2 ), where 
a = 0.5,0.421,0.241,0.126 for Z = 0,0.25,0.55,0.83, respectively. Maxima and minima 
in the differential conductance of figure 2 are associated with various crossings in figure 
3. 
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Consider for example the Z = 0.55 and Z = 0.83 results, where the first 
maximum in G~^dl / d{fii — ^2) corresponds to the crossing fj,± — fx — A s . For these struc- 
tures, provided 11-/12 < A*-, excitations from the right reservoir (2) are almost completely 
Andreev reflected at x = L and therefore the conductance is dominated by scattering of 
electrons from the left reservoir at x — 0. In this limit, it is of interest to examine the 
quantity Gj f 1 dl/d(fj,i — //), which in the absence of quasi-particle transmission, is equiv- 
alent to the left boundary conductance, examined by BTK. The solid lines of figure 4 
show self-consistent results for this quantity. For comparison, the dashed lines show non 
self-consistent results (which are essentially those of BTK [17]) obtained by insisting that 
\i = fj, 2 and for < x < L, A(x) = A For Z = 0.83 the dashed and solid curves of 
figure 4 are in good agreement, reflecting the fact that for large Z the current is small 
and therefore A(x) is not significantly modified. For the smaller barrier strength, the 
self-consistent conductance differs significantly from the BTK curve. In the presence of 
quasi-particle transmission [18,20], the resistance of a N-S-N structure does not reduce 
to the sum of two boundary resistances. Consequently even for Z = 0.83, the N-S-N dif- 
ferential conductance of figure 2 shows extra structure which is absent from a boundary 
conductance calculation. For example the second peak of the Z = 0.83 results of figure 2 
corresponds to the crossing /i — /i 2 = A s — v s pf, while the minimum between these peaks 
corresponds to the crossing ^ — fi = A s + v s pp. For Z = 0.55 the peaks at /ii — /i and 
/i — /i 2 are no longer separated, but again a minimum occurs at /ii — fj, — A s + v s pf- For 
this value of Z, a maximum occurs at /ii — /U 2 ~ 3A S , at which the magnitude of A(L/2) 
starts to become significantly reduced by the current. 

To obtain the above results, we have presented the first self-consistent de- 
scription of a superconducting nano-structure, which incorporates quasi-particle phase 
coherence and non-equilibrium effects. The computed current-voltage characteristic of 
a single superconducting island is the result of several competing phenomena associated 
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with quasi-particle scattering from boundaries and modification of the order parameter by 
both superflow and a non-equilibrium quasi-particle distribution. This competition pro- 
duces significant new structure, particularly at energies above A , which is not contained 
within a non self-consistent description. This structure is observable experimentally [19] 
but in earlier discussions [15] has been dismissed as a distortion due to heating. 



Note added in proof 

A similar calculation has been carried out recently by F.Sols and J. Sanchez- 
Canizarez, in which the superconductor is treated as an incoherent reservoir [21]. Where 
agreement is expected, their results are comparable with those reported here. 
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Figure captions 



Figure 1. Results for the case Z = 0. The upper figure shows the form of the at- 
tractive interaction V(x) (in units of // ) used in all the calculations of this paper. 
The middle and lower graphs show the self-consistent forms the magnitude and phase 
of A(x) = \A(x)\expi<f>(x). The results with the largest values of |A(x)| and constant 
phase correspond to /ii — fj, 2 — 0. In order of decreasing |A(x)| and increasing phase 
gradient, remaining results correspond to ^1 — ^2 = 0.006 and //1 — //2 = 0.01, respectively. 

Figure 2. Self-consistent results for the differential conductance of a superconducting 
dots, with a delta function barrier of strength Z, located at x — 0. 

Figure 3. For each value of Z, the upper (lower) dashed line shows self-consistent 
results for the voltages //+ = A s +Pfv s (/i- = A s —Pfv s ), while the thin solid line shows 
A s = |A(L/2)|. The upper (lower) thick solid line shows self-consistent values of /ii — /i 
(/i — /i 2 ). For Z — 0, the latter are equal. 

Figure 4. The solid lines show self-consistent results for Gj^dl /d([ii — //). The dashed 
lines show non self-consistent results obtained by setting A(x) = A s for < x < L and 
zero elsewhere. 
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